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Abstract 



Large-space no-core shell model calculations have been performed for 3 H, 4 He, 
5 He, 6 Li, and 6 He, using a starting-energy-independent two-body effective in- 
teraction derived by application of the Lee-Suzuki similarity transformation. 
This transformation can be performed by direct calculation or by different 
iteration procedures, which are described. A possible way of reducing the 
auxiliary potential influence on the two-body effective interaction has also 
been introduced. The many-body effects have been partially taken into ac- 
count by employing the recently introduced multi-valued effective interaction 
approach. Dependence of the 5 He energy levels on the harmonic-oscillator 
frequency as well as on the size of the model space has been studied. The 
Reid 93 nucleon-nucleon potential has been used in the study, but results have 
also been obtained using the Nijmegen II potential for comparison. 



I. INTRODUCTION 

Large-basis no-core shell-model calculations have recently been performed [l|-|8| . In these 
calculations all nucleons are active, which simplifies the effective interaction as no hole states 
are present. In the approach taken, the effective interaction is determined for a system of 
two nucleons only and subsequently used in the many-particle calculations. To take into 
account a part of the many-body effects a multi-valued effective interaction approach was 
introduced and applied in the no-core shell-model calculations || and also tested in a model 
calculation |J. 

In these shell-model calculations different approaches have been taken in deriving the two- 
nucleon effective interaction. All the methods employed, however, relied on the reference 
G-matrix method, introduced in Ref. (|10f| , which leads to two-body matrix elements of 



a starting-energy-dependent G-matrix. To get rid of this unwanted dependence either a 
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suitable parametrization was chosen f§-§],|[ or folded diagram effects were taken into account 
by calculating the derivatives of the G-matrix in an approximate way 



In the present paper we apply the Lee-Suzuki similarity-transformation approach [ITJ to 
derive the two-body effective interaction. We try to avoid unnecessary approximations in 
performing the calculations. The harmonic-oscillator insertions are kept; consequently, the 
effective interaction is A-dependent. Also the hermitization of the effective interaction, which 
is, in general, non-hermitian, is done without approximations by a similarity transformation. 

We consider three possible ways of deriving the effective interaction. The first two are the 
standard iterative procedures, starting with the G-matrix calculation pD| . We study the Lee- 



Suzuki vertex renormalization iteration, which makes use of the G-matrix derivatives. Unlike 
most of the previous applications of this approach, we calculate these derivatives exactly 
employing the reference G-matrix derivatives. The other iterative procedure, usually refered 
to as Krenciglowa-Kuo technique, is carried out by diagonalizing the non-hermitian effective 
hamiltonian in subsequent iterations ||12|| . This method has so far been applied only to model 
calculations. We also discuss its recent generalization fL3f . Finally, we directly construct 
the effective interaction without calculating the G-matrix by employing the solutions of the 
two-body problem. When criteria necessary for convergence are the same for both iteration 
procedures, the resulting effective interactions obtained in all three methods are identical. 
The last method, however, has two advantages; first, its simplicity, and second, the fact 
that it utilizes an explicit construction of the transformation operator. This transformation 
operator may then be used for the calculation of other effective operators. 

To take partially into account the many-body effects neglected when using only a two- 
body effective interaction, we employ the recently introduced multi-vauled effective interac- 
tion approach ||. 

It was observed earlier J7| that the two-body effective interaction derived by the vertex 
renormalization method is too attractive and leads to overbinding of the many-body system. 
We discuss this problem and believe that it is largely due to the uncompensated Q space 
part of the auxilliary harmonic-oscillator potential. We discuss possible treatment of this 
problem by modification of that part of the auxilliary potential. 

In section [TT] we discuss the shell model hamiltonian with the bound center-of-mass as 
well as the two-particle hamiltonian and the methods used to derive the starting-energy- 
independent effective interaction. Results of the calculations for 3 H, 4 He, 5 He, 6 Li, and 6 He 
are presented in section |TJ. In particular, we discuss the harmonic-oscillator frequency and 
the model-space-size dependences of the 5 He states. Conclusions are given in section [TV]. 

II. SHELL MODEL HAMILTONIAN AND STARTING-ENERGY 
INDEPENDENT EFFECTIVE INTERACTION 

In most shell model studies the one- and two-body hamiltonian for the A-nucleon system, 

i.e., 

1=1 z " i i<j 

where m is the nucleon mass and Vij the nucleon-nucleon interaction, is modified by adding 
the center-of-mass harmonic oscillator potential |v4mf2 2 i? 2 , R = \ j2f=i^i- This potential 
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does not influence intrinsic properties of the many-body system. It provides, however, a 
mean field felt by each nucleon and allows us to work with a convenient harmonic oscillator 
basis. For an alternative manipulation of the center-of-mass terms see, e.g., Ref. M. The 
modified hamiltonian, depending on the harmonic oscillator frequency Q, can be written as 

A 

1=1 

which is the same as Eq. (4) in Ref. [|lj . Shell-model calculations are carried out in a model 
space defined by a projector P. In the present work we will always use a complete NhQ 
model space. The complementary space to the model space is defined by the projector 
Q = 1 — P. Consequently, for the P-space part of the shell-model hamiltonian we get 

P. (3) 

off 

The effective interaction appearing in Eq.(^) is, in general, an A-body interaction, and, if it is 
determined without any approximations, the model-space hamiltonian provides an identical 
description of a subset of states as the full-space hamiltonian (||). The intrinsic properties 
of the many-body system still do not depend on Q. From among the eigenstates of the 
hamiltonian (|3|), it is necessary to choose only those corresponding to the same center-of- 
mass energy. This can be achieved by projecting the center-of-mass eigenstates with energies 
greater than ^hQ upwards in the energy spectrum. The shell-model hamiltonian, which does 
this, takes the form 

H%= E P 

i<j=i 

+ pp(H? m -hn)p, (4) 

where (3 is a sufficiently large positive parameter and = £f!i;+ \AmVt 2 R 2 , P cm = Y^t=i Pi- 
In a complete NhQ model space the removal of the spurious center-of-mass motion is exact. 
When going from Eq.(|3|) to (|J), we added (f3 — \)PH^ m P and subtracted /3§MXP, which has 
no effect on the intrinsic spectrum of states with the lowest center-of-mass configuration. 

The effective interaction should be determined from H (0). Calculation of the exact A- 
body effective interaction is, however, as difficult as finding the full space solution. Usually, 
the effective interaction is approximated by a two-body effective interaction determined from 
a two-nucleon problem. The relevant two-nucleon hamiltonian obtained from (0) is then 

fi? = Hg + V? = + + rl) + V(n - f 2 ) - ^(rl - f 2 ) 2 . (5) 

This can be transformed to two-nucleon relative and center-of-mass parts by introducing the 
coordinates q = — P2), P2cm = Pi + P2, R2cm = + f 2 ) and r = fx — f 2 yielding 

«? = % + \M^» + % + 4^Vi* + V(r, , (6) 
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with M = 2m and \x = \m. While the center-of-mass part has the solution Ej^ c = (2A/"+£ + 
|)M] and the eigenvectors \Af£), the relative-coordinate part can be solved as a differential 
equation or, alternatively, can be diagonalized in a sufficiently large harmonic oscillator 
basis. The latter possibility is, obviously, not applicable for hard-core potentials. 

The starting-energy-dependent effective interaction or G-matrix corresponding to a two- 
nucleon model space defined by the projector P 2 can be written as 

G(e) = V? + V?Q* e _Q,H2 Q Q*V? > ( 7 ) 

where Q 2 = 1 — P2 and V 2 is the interaction given by the last two terms on the rhs of 
Eq.(||). The G-matrix (|7|) can be constructed from the solutions of the Schrodinger equation 



with the hamiltonian d^), by using the reference matrix method ||10|| . The G-matrix can be 
expressed as 

G(e) = A(e)- 1 G R (e) , (8a) 

Gn(e) = (Hg - a) + (fl£ - e) £ j^« 2 - e) , (8b) 

k £ ~ ^ k 

A(e) = 1 + G K {e)P 2 —^ , (8c) 

where Hq 2 is given by the first two terms on the rhs of Eq.([5[) and E^, \k) are the eigenvalues 
and eigenvectors of H 2 (|^), respectively. 

To obtain a starting-energy-independent effective interaction, one has to take into ac- 
count the folded diagrams or, equivalently, to construct a similarity transformation that 
guarantees decoupling between the model space P and the Q-space. We employ the Lee- 
Suzuki [11] similarity transformation method, which gives the effective interaction in the 
form 

P2V 2cS P 2 = P 2 V 2 n P 2 + P 2 V 2 n Q 2 uP 2 , (9) 
with lo satisfying the equations u> = Q 2 tuP 2 and 

u = & nunn Q2V2P2 ~ Q2 n \ nn Q2toP 2 (H? + H%Q 2 u - s)P 2 . (10) 

e — LJ 2 U 2 LJ 2 e — <^ 2 n 2 LJ 2 

In this degenerate-model-space formulation two iterative solutions of the equation (|i0|) exist 
and lead to different expressions for the effective interaction. The first one, the Krenciglowa- 
Kuo (KK) iteration procedure [0, gives the for the n-th iteration formula 

V 2eS , n = £ P 2 G(e + E l}n ^)P 2 \l n ^)(l n ^\P 2 . (11) 
1 

In ( [TT] ) the states |/ n -i) are the right eigenvectors of Hq 2 + V 2 eS,n-i — £ belonging to the 
eigenvalue Ei >n -i. The tilda states are the biorthogonal eigenvectors. When this procedure 
converges, it does so to the states which have the largest overlap with the model-space 
states. The resulting V 2e s © is independent of the starting energy e. This method was 



recently generalized to a non-degenerate model space | 13| . The difference is in the starting 
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iteration. When using (pTTj) , the starting iteration is G(e), while the non-degenerate model 
space version |H| starts with 



2a 



(12) 



Here the e a are the unperturbed energies, in our case the two-nucleon harmonic-oscillator 
energies, and P 2a is the projector on the two-nucleon state \a). Unlike in the original paper 
T3| , we introduce a shift A, as the G-matrix cannot be evaluated for the energy e a , using 



the reference matrix method (g). Note that fll2|) is the interaction used in the previous 
no-core shell-model studies with A treated as a free parameter. 

One can also use the alternative procedure, usually called the vertex-renormalization 
approach, which can be obtained from Eq.([TT]) |TTJ. The resulting iteration sequence for the 
effective interaction can be written as 



Hq 2 + V 2e S,n — £ 



-E 



2 — 

n-1 



P2G {1 \e)P 2 
1 



m=2 



ml 



P 2 G {m \e)P 2 (Hg + y 2eff ,„- m+ i - e)P 2 (Hg + V 2efr , n _ 



m+2 



-1 



+ V 2eS , n -i - e)P 2 P 2 + G(e) - e) P 2 



-e)P 2 



(13) 



where V2eff,i — P2G(s)P 2 . This method, which relies on the G-matrix derivatives G^ m \ was 
applied in shell-model calculations for the first time by Poppelier and Brussard In 
that and most other applications different numerical approximations were used to evaluate 
the derivatives PJ7|JI^]. In fact, these derivatives can be calculated exactly, because the 

) can be differentiated to any 



reference matrix Gr(e) (|8b ) as well as the operator A(e) ( 
order analytically, as can be seen when they are expressed in the two-nucleon harmonic- 
oscillator basis. For example, the n-th derivative, n > 0, of Gr is 



Gr2 7 (^) 



'-l) n n\ 



+ (e - £ a ) 
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(14) 



and similarly for A(e). Then the n-th G-matrix derivative can be expressed as 

n | 

G (n) = ^ _ ^__ ( jA -l)(n- ro ) G rM > 



m=0 



n 



m M 



(15) 



where the derivatives of the inverse matrix can be obtained from 



-l\(n) 



n—l 

E 



nl 



_^-l^(n-m)^-l^(m) 



(16) 



m=0 (n-m)!m! 

An alternative way of calculating G-matrix derivatives analytically was proposed in Ref. 



TBI . When the convergence of (]T~3|) is achieved, the effective interaction does not depend on 



the starting energy e. The states reproduced in the model space are those lying closest to e. 
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As our calculations start with exact solutions of the hamiltonian ([|), we are, in fact, 
in a position to construct the operator oj and, hence, the effective interaction directly from 
these solutions. Let us denote the two-nucleon harmonic-oscillator states, which form the 
model space, as \cep), and those which belong to the Q-space, as \o>q). Then the Q-space 
components of the eigenvector \k) of the hamiltonian (|5|) can be expressed as a combination 
of the P-space components with the help of the operator uj 



(ot Q \k) = ^2(a Q \u\a P )(a P \k) 



(17) 



(Xp 



If the dimension of the model space is dp, we may choose a set K, of dp eigenevectors, 
for which the relation ( |T7j ) will be satisfied. Under the condition that the dp x dp matrix 
(otp\k) for \k) G /C is invertible, the operator uj can be determined from flT7|). Note that in 
the present application the eigenvectors \k) are direct products of the center-of-mass and 
the relative coordinate eigenvectors. The condition of invertibility is not satisfied for an 
arbitrary choice of the eigenvector set /C. Once the operator uj is determined the effective 
hamiltonian can be constructed as follows from Eq.(RI) 



(lp\H 2eS \a P ) 



E 



(j P \k)E k (k\a P ) + ^2^p\k)E k (k\aQ){a Q \uj\ap) 



(18) 



It should be noted that P 2 \k) = J2 ap \oip){otp\k) for \k) G fC is a right eigenvector of ( |18l) 
with the eigenvalue E k . 

In the case when the iteration conditions are the same for both methods fllTD and (|13|). 
and the set /C, for which Eq.([T7|) is fulfilled, is chosen accordingly, all three methods lead to 



the identical effective hamiltonian. This hamiltonian, when diagonalized in a model space 
basis, reproduces exactly the set JC of dp eigenvalues E k . Note that the effective hamiltonian 
is, in general, non-hermitian, or more accurately quasi-hermitian. It can be hermitized by a 
similarity transformation. When the direct method ( TT7|JTE| ) is used the similarity transfor- 
mation is determined from the metric operator P 2 (l + co*u))P2. The hermitian hamiltonian 
is then given by [T(J 



H 



2eff 



P 2 (l+Ju)P 2 



H- 



2cff 



P 2 (l + Juo)P 2 



(19) 



When the iteration methods (|TTD or ( |13|) are employed, the u operator is not determined. In 



previous applications PJ7|fl^|, the hermitization was usually done by averaging the conjugate 
matrix elements. However, also in this case a similarity transformation can be constructed, 
which hermitizes the effective hamiltonian. As shown in Ref. [17], if S is a matrix, which 
diagonalizes H 2eS , SH^S' 1 = D, then the metric operator can be expressed as S^S. 
Consequently, it follows that the hermitized H 2e s is obtained from 



2eff 



S^S 



H- 



2cfT 



s*s 



(20) 



Finally, the two-body effective interaction used in the shell-model calculation is deter- 



mined from the two-nucleon effective hamiltonian as V 2e fi = H 2c g — H, 



02- 
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III. APPLICATION TO LIGHT NUCLEI 



In this section we apply the methods for calculating the two-body effective interaction 
outlined in section |Tj and, with the obtained interactions, we perform no-core shell-model 
calculations for nuclei with A = 3 — 6. We use a complete NhQ model space with, e.g., N = 8 
for the positive-parity states of 4 He. This means that 9 major harmonic-oscillator shells may 
be occupied in this case. The two-nucleon model space is defined in our calculations by iV max , 
e.g., for an 8hQ calculation for 4 He A max = 8. The restriction of the harmonic-oscillator shell 
occupation is given by N\ < N max , N 2 < N max , (Ni + N 2 ) < iV max . The same conditions 
hold for the relative 2n + l, center-of-mass 27V" + £, and 2n + l + 2J\f+jC = Nx + N%, quantum 
numbers. 

In the present calculations we use the Reid 93 nucleon-nucleon potential [|18|| . For 4 He 
we make a comparison with the Nijmegen II potential [Tj| with corrected l Pi wave [HJ. We 
work in the isospin formalism and do not include the Coulomb interaction. The np channel 
interactions of Reid 93 and Nijmegen II are used. 

First, let us compare the different methods for calculating the starting-energy- 
independent effective interaction. The KK method (|TT| ) works well for iV max < 4. Then 
the convergence is achieved usually after less than 15 iterations and the model space eigen- 
values differ from the full space eigenvalues, which are of the order of 10 1 — 10 2 MeV, by not 
more than 10~ 4 MeV. Note that the reference G-matrix (|8b|) becomes singular for e = E^. 
Consequently, exact eigenvalues cannot be reproduced when the reference G-matrix method 
is applied for calculating the G-matrix. However, this singularity causes no problem when 
the iteration is stopped after achieving the above mentioned precision. When the starting 
iteration ([12]) is used instead of G(e), usually a few iteration steps are saved. As to the 
starting energy, e = is a possible choice. In most calculations we used negative values for 
e. For A a non-zero value must be chosen, we used typically about -5 MeV. The resulting 
effective interaction is not dependent on these choices, and they also do not affect the num- 
ber of iterations significantly. For iV max > 4 divergence in some channels can be encountered 
and, moreover, the calculation becomes time consuming. 

The application of the vertex-renormalization method ( ]TB| ) requires the G-matrix deriva- 
tives. The matrix elements of the G-matrix derivatives decrease rapidly but, on the other 
hand, they are multiplied by effective-interaction matrices of increasing powers. Conse- 
quently, the overall convergence is rather slow for larger model spaces. Also the rate of 
convergence differs for different states. When the starting energy is chosen below the ground 
state energy, the fastest convergence is achieved for the lowest states, which have the largest 
overlap with the model space. This is the same observation as found in the model calcu- 
lations pup . The rate of convergence is very sensitive to the choice of the starting energy 
e. The closer to the ground state, the faster the convergence. On the other hand, serious 
numerical problems occur when a higher number of iterations is required. These problems 
are related to the above mentioned fact that we multiply very small numbers by very big 
numbers when calculating higher iterations. To achieve the same accuracy as with the previ- 
ous method, often more than 30 iterations are needed, and to curb the numerical difficulties 
real* 16 precision is neccessary. These facts make this method rather impractical. It is also 
difficult to apply this method for model spaces with A max > 4. 

Let us point out that new iteration methods were suggested recently, which combine 
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some features of both methods studied here These techniques are more involved from 



the computational point of view, and we did not try to use them. They may, however, 
remedy some of the difficulties found in our applications. 

The simplest and the most straightforward way to calculate the two-body effective in- 
teraction is accomplished by solving the equation ( |T7| ) and constructing the effective hamil- 
tonian according to Eq . (|18D . The easiest way to perform the calculations is in the center- 
of-mass and relative coordinate basis and afterwards to do the transformation to the two- 
nucleon harmonic-oscillator basis. Note that uj is diagonal in the center-of-mass quantum 
numbers Af, £ as well as in S, j, J, T. Consequently, the sum in Eq.flT7|) goes only over rip, lp 
for the basis states classified by \nlSj,M£,JT). An important point is the right choice 
of the eigenstates in the set K. For each S,j,Af,£, J,T we must choose as many relative 
coordinate eigenstates as the allowed number of rip, lp combinations. These are determined 
from the condition 2Af + £ + 2rip + lp < N max . Since the harmonic-oscillator basis is infinite, 
we make a truncation in the Q-space by keeping only the states with tiq < 150. Note that 
I — j ± 1 for the coupled channels and I = j for the uncoupled channels. This method can 
be easily applied to calculate the effective interaction up to N max = 8, as required in the 
present shell-model calculations. 

To take partially into account the many-body effects neglected when using only a two- 
body effective interaction, we employ the recently introduced multi-valued effective interac- 
tion approach ||. In this approach different effective interactions are used for different hfl 
excitations. The effective interactions then carry an additional index indicating the sum of 
the oscillator quanta for the spectators, N sps , defined by 

M = N - N = N' - N (21) 

iv sps J, sum ly a iv sum 7 ' \ ) 

where N snm and N' Bwa are the total oscillator quanta in the initial and final many-body states, 
respectively, and N a and are the total oscillator quanta in the initial and final two-nucleon 
states \a) and I7), respectively. Different sets of the effective interaction are determined for 
different model spaces characterized by N sps and defined by projection operators 

P2{N eps ) = l-Q 2 (N sps ) . (22b) 

This multi-valued effective-interaction approach is superior to the traditional single-valued 
effective interaction, as confirmed also in a model calculation ||. 

The shell-model diagonalization is performed by using the Many-Fermion-Dynamics 
Shell-Model Code fl2~2"f . We present the results for A = 3 — 6 nuclei in Tables It 



has been observed earlier and it is apparent in the present calculations as well that, 
when the effective interaction is derived using the Lee-Suzuki method, the interaction be- 
comes too strong and leads to overbinding of nuclei. The problem is not in the calculational 
method but rather in the fact that the original two-nucleon hamiltonian (^) is flawed if 
P2 ^ 1. The difficulty is that while the relative-coordinate harmonic-oscillator auxiliary 
potential is exactly cancelled in the equation @, it is not fully cancelled when the effective 
interaction is derived from the two-nucleon hamiltonian (|5|) . This can be seen, for example, 
when the nucleon-nucleon interaction is switched off. Then from Eq.(|]) we obtain an ef- 
fective interaction derived from the relative-coordinate harmonic-oscillator potential, which 



S 



will be different from the model-space part of the harmonic-oscillator potential, appearing 
in Eq.(|J). Clearly, the Q 2 -space part of the relative-coordinate harmonic-oscillator potential 
is responsible for this effect. In order to reduce this spurious effect of the auxiliary potential 
on the effective interaction, we introduce an ad hoc modification of the relative-coordinate 
two-nucleon Q-space part of the auxiliary potential as follows 

= £ + P 2 ^Vr-P 2 + P 2 k Q ^^Q 2 

+Q 2 k Q ^^^P 2 + Q 2 k Q ^^^Q 2 



+ (l-k Q )hn ] j^j^(N max + 2)Q 2 + V(r) . (23) 

Here we have introduced a constant kq < 1; for kq = 1 we get the original hamiltonian 
appearing in Eq.(^). Moreover, for A = 2 the harmonic-oscillator dependence vanishes and 
for P — > 1 the relative-coordinate part of the hamiltonian (|6]) is recovered. In Eq. (|23|) the 
Q-space harmonic-oscillator potential is made more shallow and is shifted. The shift is such 
that the P-space and the Q-space potentials are equal for r 1 = |((r e )jv max + i^ 2 ) N max +i) , 
with this mean value determined for the eigenstates of the relative-coordinate harmonic- 
oscillator hamiltonian appearing in Eq. (|^). In this modified hamiltonian (p3|), a quasi- 
particle scattered into the Q-space feels a weaker auxiliary potential. Note that an alternative 
method by Krenciglowa et al. [^8| for deriving the G-matrix, as opposed to that we employ 
here uses plane- wave- type intermediate Q-space states, that is, no auxiliary potential 



at all. For a recent review of this approach see Ref. [2~9~|. Our present modification, in 
fact, brings those two methods closer together. We solve the Schrodinger equation with the 
hamiltonian ( PB"|) by diagonalization in a harmonic-oscillator basis characterized by b = 
with the radial quantum number n = . . . 150. The error caused by this truncation of 
the harmonic-oscillator basis can be estimated for kq = 1, when the system can be solved 
as a differential equation. We found that the low-lying eigenvalues obtained in the two 
calculations do not differ by more than 10 -3 MeV and in most cases by much less. The 
lowest eigenvalues are typically of the order of 10 1 MeV. The diagonalization cannot be 
used for hard-core nucleon-nucleon potentials, but for the soft-core Reid 93 and Nijmegen 
II potentials the interaction matrix elements can be evaluated straightforwardly. Note that 
the hamiltonian (|23|) depends on iV max and the projection operators. When the multi- 
valued effective interaction is calculated, solutions are found only for projectors and iV max 
corresponding to N sps = 0. The value iV max = 8 is used in most calculations. Using these 
solutions, all required sets of effective interactions are constructed. 



The question arises if the modification we introduced in Eq. fl23|) may not cause center- 
of-mass spurious contamination of the physical states. We note that the modification is done 
only in the Q-space part of the auxiliary potential and, most importantly, with regard to 
the two-nucleon relative coordinate. In conjuction with this, our choice of the two-nucleon 
model space as well as the complete NhQ model space for the many-nucleon states prevents 
the center-of-mass contamination of physical states. We have tested numerically for possible 
spurious center-of-mass contaminations by varying the parameter ft, introduced in Eq.(|j), 
in calculations for several systems and found that the physical states remain unchanged for 
different choices of (3, including (3 — 0, even if kg ^ 1. 
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We note that the free (or bare) values of the nucleon charges are used for calculating 
mean values of different operators. 

In Table | we present results of the 8hQ calculation for 3 H with Ml = 19.2 MeV as 
suggested by the formula |3(J (in units of MeV): 

hQ = 45A~^ - 25A'I . (24) 

We observe that the calculation, using the effective interaction derived from the hamiltonian 
or equivalently Icq — 1 in Eq. (p3|), overbinds the system in comparison with both the 
experimental value (8.482 MeV) and the exact result for the Reid 93 calculation (7.63 MeV) 
|3jfl . Our aim should be to reproduce the result of the exact calculation. To achieve that 



we varied the parameter kQ in Eq. (|23|). The value fcq = 0.6 gives reasonable agreement. 
We keep this value for all other calculations, so that a meaningful comparison of binding 
energies can be obtained. 

In Table |I| the calculated results for 4 He are presented, where we have used 8hQ for the 
positive-parity states and 1MI for the negative parity states. The value Ml = 18.4 MeV, 
obtained from fl2~4|) is employed. The results for calculations with kQ — 1 and fcg = 0.6, 
using the Reid 93 potential are shown, as well as calculations with fcq = 0.6, using the 
Nijmegen II potential. We observe that both potentials give almost identical results. The 
effective interaction derived using the unmodified relative-coordinate two-body hamiltonian 
overbinds the system. On the other hand, the calculation with the modified hamiltonian 
gives very reasonable agreement with the experimental results. Unlike the experimental 
spectrum we still get the 2hQ + state above the lhQ 0~; however, the discrepancy is 
reduced considerably in comparison with previous calculations || and most of the states 
compare better with the experiment data than in the recent calculation with multi-valued 
effective interaction ||. 

The calculation for 5 He was performed for several values of hfl and different model spaces 
in order to study the dependence of different states on these quantities. In Table |TTT| we show 
the 1MI (n = +) and QMl (n = — ) results, respectively, using Ml = 17.8 MeV obtained from 
(E4j). The controversy regarding the shell-model calculations for this nucleus has to do with 



the nature of the excited states 0]. In the standard shell- model formulation, also employed 
here, the center-of-mass of the nucleus is bound in a harmonic-oscillator potential, so all 
states are bound. However, they do not necessarily belong to the internal excitations of the 
studied nucleus. They may as well correspond to, e.g., a two-cluster configuration. One 
would expect that such states are more sensitive to the variation of the model space size and 



the harmonic-oscillator binding potential |3"3 ]. In Fig. |Tj we present the hfl dependence of 
the 5 He states calculated in the 7MI (n = +) and 6hQ (ir = — ) model spaces, respectively, 
by using the effective interaction obtained for /cq = 0.6 and iV max = 8 in Eq.(p3f). Note 
that only three experimental states are known in this nucleus. Also the excitation energy 
of the | state is determined with a large error. For comparison, in Fig. ^| results for the 
5hQ (n = +) and 4hQ (tt = — ) model spaces, respectively, are presented. The effective 
interaction for this calculation has been derived from Eq . (^3|) with the same /cq = 0.6 but 
different iV max = 6. We observe large sensitivity for the higher excited states on changes in 
Ml with the exception of | 2 . This state, dominated by the (0s) 3 (0p) 2 configuration, is a 

good candidate for the experimental | + state. Note also the significant shifts of the excited 
states, again with the exception of | 2 and ^ , downwards in energy, when the model space 
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is enlarged. Even if we cannot draw a definitive conclusion, these facts indicate that excited 
states obtained in the calculation, but unobserved, may be states, which do not correspond 
to 5 He internal excitations. 



The calculated results for Li are presented in Table [TV]. A 6M7 calculation was per- 
formed for the positive-parity states using Ml = 17.2 MeV obtained from (|24|). Again the 
effective interaction, derived from the hamiltonian @, is too strong and overbinds the sys- 
tem. The calculation with the modified hamiltonian (E3T) gives very reasonable agreement 
with experiment for all properties shown. In addition to the experimental states presented 
in the Table [TV], another J 71 = 3 + , T = state with the excitation energy 15.8 Mev is known. 
The lowest J n = 3 + ,T = 2MI state obtained in our kg = 0.6 calculation has an energy 
of 13.840 MeV, with a main configuration of (0s) 4 (0p) 1 (lp0/) 1 . Another candidate for this 
experimental state may be the calculated state with an excitation energy of 18.378 MeV, 
dominated by the configurations (Os) 4 (lsOof) 2 and (0s) 3 (0p) 2 (ls0<i) 1 . Likely, the latter state 
may be more stable against model space and hQ variations. We also performed a 5hQ cal- 
culation for the negative-parity states of 6 Li. The lowest calculated state is J 71 " = 2~, T = 
with an excitation energy 9.135 MeV, followed by J n = 1~,T = (£^=9.335 MeV) and 
J n = 0~,T = (-£^=11.372 MeV). Such states have not been observed experimentaly. 

In Table [V] we present the 6 He characteristics. As isospin symmetry is not broken in the 
present calculations, the energies are the same as those for 6 Li with T = 1. In addition, 
the proton and neutron radii are evaluated. Reasonable agreement with the experimantal 
values is found. 

Contrary to the previous no-core shell- model calculations ||[8|, the effective interaction 
used in the present calculations does not contain a parameter A adjusted for each nucleus 
in order to get a reasonable binding energy. Besides the model-space size, the present 
calculations depend only on the harmonic-oscillator frequency hfl and on the parameter kQ, 
appearing in the relative coordinate two-nucleon hamiltonian (|23"D , which distinguishes the 
structure of the P-space from that of the Q-space. For the results presented in Tables UW\ we 
chose Ml according to formula (|24"D and kQ was fixed to reproduce the result of the exact 3 H 
calculation. Consequently, our calculations contain no variable parameters and meaningful 
comparisons can be made between our calculated binding energies and the quantities derived 
from them and the experimental values as well as the results of other calculations. 

In this regard, we note the recent calculation of valence energies of 6 He and 6 Li and the 
neutron separation energy of 5 He, using the two-frequency shell-model approach p4 |. For 



example, our calculation with kq = 0.6 and the Reid 93 potential gives for the neutron 
separation energy of 5 He, E sp = i? B ( 5 He) — i?B( 4 He) = —1.303 MeV to be compared with 
the experimental value of —0.894 MeV. Similarly, the valence energy of 6 He is found to 
be £ va i( 6 He) = -[£ B ( 6 He) + E B ( 4 Re) - 2E B ( 5 B.e)} = -3.151 MeV, compared with the 
experimental -2.761 MeV. Furthermore, £7 vai ( 6 Li) = -[E B ( G U) + £ B ( 4 He) - E B ( 5 Re) - 
E B ( 5 U)] = -6.845 MeV, where we used E B ( 5 Li) = E B ( 5 Re) = 26.105 MeV from Table 0, 
as the Coulomb contributions, not taken into account in the calculations, more or less cancel. 
Here the experimental value is —6.559 MeV. We observe good agreement of these quantities 
with the experimental values, in fact, better than that obtained in Ref. |j34 |. Apparently, 



the reason for these improved results is the size of the model space, which is much larger in 
our calculations. 

In Ref. |32[] concerns were raised about the stability of 6 Li to the a + d threshold in 
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previous no-core calculations From the present results we deduce that 6 Li is bound 
against this threshold by 1.314 MeV, compared with the experimental value of 1.474 MeV. 
Note that the deuteron is treated exactly in our formalism, Es{d) = 2.2246 MeV. To arrive 
at the present number, we used the Coulomb energy contributions to binding energy Ec{Z = 



2) = —0.76 MeV and Ec(Z = 3) = —1.46 MeV. Another issue raised in Ref. [32] was sign 
of the quadrupole moment of 6 Li. As in the previous no-core calculations we also get 
the sign correctly in our Icq = 0.6 calculation. Clearly, this is a consequence of a reasonable 
effective interaction and a large enough model space. 



IV. CONCLUSIONS 

In the present paper we have presented different ways of constructing a starting-energy- 
independent two-body effective interaction. In particular, we studied the Lee-Suzuki sim- 
ilarity transformation method and compared two different iteration schemes for obtaining 
a solution, the vertex renormalization method and the Krenciglowa-Kuo method. When 
applying the vertex renormalization method, we obtained the derivatives of the G-matrix 
exactly by using the derivatives of the reference G-matrix. A third approach involved a 
direct calculation of the transformation operator u. When the convergence conditions of 
the iteration methods are satisfied, all three approaches lead to the same two-body effective 
interaction. Our conclusion is that, for large model spaces that we employ, the only viable 
option is the direct calculation method. For large model spaces (iV max > 4), the iteration 
methods either fail to converge or the calculations are time consuming or both. 

Unlike in most previous applications, where the non-hermitian effective interaction was 
hermitized by averaging the conjugate matrix elements, we hermitize the effective interaction 
exactly using a similarity transformation. It was demonstrated in Ref. pS that the averaging 



is a good approximation for the hermitian sd and pf effective interactions. We observe here, 
however, that for large model spaces, like those we employ, the non-hermiticity could be 
significant. It is certainly preferable to work with the exactly hermitized interaction. 

Employing the derived effective interaction, we have performed shell-model calculations, 
in which up to 9 major harmonic-oscillator shells may be occupied. In order to take into 
account part of the many-body effects, we utilized the new multi-valued effective interaction 
approach ||. The results for no-core, full NhQ calculations are given for nuclei with A = 
3 — 6. The effective interaction was derived from the Reid 93 nucleon-nucleon potential, 
but calculations were also made with Nijmegen II potential for comparison. As observed 
earlier |7[], when the Lee-Suzuki method is applied for a two-nucleon system with a harmonic- 
oscillator auxiliary potential, the resulting effective interaction is too strong and leads to 
overbinding of the many-body system. This problem is caused by incomplete cancellation 
of the relative coordinate part of the auxiliary potential. To mend this flaw, we introduced 
a modification of the Q-space part of the relative-coordinate two-nucleon hamiltonian, from 
which the effective interaction is calculated. In effect this weakens the auxiliary potential in 
the Q-space. 

In the present calculations, besides the model space size, the only free parameters are the 
harmonic-oscillator frequency hfl and the parameter modifing the two-nucleon hamiltonian, 
as discussed above. The latter parameter was fixed from the 3 H binding energy calculation 
(fitted to the result of an exact 3 H calculation) and HQ was taken from the phenomenological 
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formula (p4|) . Hence, unlike previous no-core calculations, we are able to compare, for nuclei 
other than 3 H, quantities derived from binding energies, such as valence energies, with 
experiment. 

For most calculated characteristics we found good agreement with the experimental 
values. In agreement with the previous observation [0, we found that the Reid 93 and 
Nijmegen II potentials give very similar results. The question of low-lying positive-parity 
states in 5 He was also investigated. We observed that the calculated low-lying positive-parity 
states have not converged, regarding changes in the model-space size and variations in HQ. 
This may indicate that they do not correspond to internal excitations of 5 He. However, we 
are not in a position to make a conclusive statement concerning this issue. 

Because we derive the transformation operator u explicitly in our calculations, we are 
able to use it for calculating any effective operator, employing the approaches discussed in 
Refs. |36|j37|| , in a similar way as in our model calculations |J. Besides calculating other 



effective operators, we also intend to extend the calculations to larger A. 
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FIGURES 



FIG. 1. Energy dependence for the 5 He states on the harmonic-oscillator energy HQ in a full 6 
and 7 HO, calculation. Only the three lowest negative-parity states are shown. The positive-parity 
states | + , | + and | + , | + have very close energies. Therefore, only single lines for these pairs of 
states are presented. The second | + state is dominated by the s 3 p 2 configuration and should 
correspond to the experimental | state. For further description of the calculation see the text. 

FIG. 2. Energy dependence for the 5 He states on the harmonic-oscillator energy h£2 in a full 4 
and 5 fitl calculation. See Fig. [I] for the details. 
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TABLES 



3 H 


k Q = l 


k Q = 0.6 


Exp 


Eb 


8.739 


7.674 


8.482 




1.540 


1.600 


1.41-1.62 




2.638 


2.634 


2.979 



TABLE I. Experimental and calculated binding energy in MeV, point proton radius in fm, 
and magnetic moment in /Un, for 3 H. The ground state has J n = The calculation 

results are presented. Different kQ choices correspond to different two-nucleon hamiltonians, from 
which the multi- valued effective interaction is calculated, as explained in the text. The Reid 93 
nucleon-nucleon potential is used, and the harmonic-oscillator parameter is taken to be HQ = 19.2 
MeV. The exact, calculated binding energy for this potential is 7.63 MeV. The free nucleon g 
factors were used. The experimental values are taken from Ref. [p3|]. 



4 He 


Reid 93 


Reid 93 


Nijmegen II 


Exp 




k Q = l 


kg = 0.6 


kg = 0.6 






31.115 


27.408 


27.499 


28.296 


Ad 


1.378 


1.434 


1.432 


1.46 


E x (0f,0) 














E x (0+,Q) 


24.009 


21.619 


21.724 


20.21 


E x (Qi,0) 


23.506 


21.290 


21.402 


21.01 


E X (2^,Q) 


25.118 


22.852 


22.948 


21.84 


E x (2^,l) 


26.544 


24.056 


24.125 


23.33 


EMI A) 


26.874 


24.263 


24.325 


23.64 


E x (li,0) 


27.763 


25.113 


25.205 


24.25 


E x (Qi,l) 


27.940 


25.247 


25.325 


25.28 


E x (l2,l) 


28.154 


25.419 


25.487 


25.95 



TABLE II. Experimental and calculated binding energy, point proton radius in fm, and exci- 
tation energies Ex^J 7 * , T) for 4 He. All energies are in MeV. The 8hQ (it = +) and 7HQ (it = — ) 
calculation results are presented. Different kg choices correspond to different two-nucleon hamil- 
tonians, from which the multi- valued effective interaction was calculated, as explained in the text. 
Both the Reid 93 and Nijmegen II nucleon-nucleon potentials are used for comparison purposes. 
The harmonic-oscillator parameter is taken to be hO, = 18.4 MeV. The experimental values are 
taken from Refs. |HjH]. 
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5 Ho 


— 1 


t\,Q — U.U 




Eb 


oU.404 


zu.lUo 


Z 1 .4UZ 




l.Ool 


i.oy ( 


M / A 
In / A 


A* 


-1.845 


1 O A A 

-1.844 


AT / A 

IN /A 


/n 


-U.4o I 


-u.4oy 


M / A 
IN / A 


^^21 ' 2/ 


A 

u 


u 


A 

u 


^(21 ' 2) 
^(21 ' 2") 


Z.OO I 


Z.UZO 


/I 4-1 
4±1 


4.U0Z 


^ 1 1 ^ 


TvT / A 


^^21 ' 2 / 

-^rr^l ' 2/ 




o.uyy 


N/ A 


9.520 


8.172 


N/A 


^(22 ' l) 


11.200 


10.413 


N/A 


^(22 ' 2) 
^ v 2 2 ' 2 / 

^ \ 2 2 ' 2 ) 


15.049 


14.025 


N/A 


20.931 


18.652 


N/A 


21.378 


19.454 


16.75 



TABLE III. Experimental and calculated binding energy, point proton radius in fm, magnetic 
moment in quadrupole moment in efm 2 and excitation energies E X (J W , T) for 5 He. All energies 
are in MeV. The 7H£l (n = +) and 6h£l (tt = — ) calculation results are presented. Different k,Q 
choices correspond to different two-nucleon hamiltonians, from which the multi- valued effective 
interaction was calculated, as explained in the text. The Reid 93 nucleon-nucleon potential is 
used, and the harmonic-oscillator parameter is taken to be HQ = 17.8 MeV. The calculated | 2 
state is associated with the experimental | + state as it is dominated by (0s) 3 (0p) 2 configuration. 
See also figures |] and |2[ The experimental values are taken from Ref. |2(| . 
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k Q = l 


kg = 0.6 


Exp 




37.532 


31.647 


31.995 




1.995 


2.097 


2.42 




0.839 


0.839 


0.822 


Q 


0.027 


-0.052 


-0.082 


^(1^,0) 











£*(3f,0) 


2.368 


2.398 


2.19 


£.(0+1) 


4.301 


3.695 


3.56 


^(2f,0) 


4.966 


4.438 


4.31 


^(2+1) 


6.878 


6.144 


5.37 


^(1^,0) 


7.232 


6.277 


5.65 


EM A) 


10.641 


9.333 


N/A 


E x {itA) 


11.442 


10.019 


N/A 


E x {lt,ti) 
EM ,1) 


11.847 


10.467 


N/A 


13.866 


12.050 


N/A 



TABLE IV. Experimental and calculated binding energy, point proton radius in fm, magnetic 
moment in /lin, quadrupole moment in efm 2 , and excitation energies E x ( J 71 ", T) for 6 Li. All energies 
are in MeV. The 6h£l calculation results are presented. Different kq choices correspond to different 
two-nucleon hamiltonians, from which the multi- valued effective interaction was calculated, as 
explained in the text. The Reid 93 nucleon-nucleon potential is used, and the harmonic-oscillator 
parameter is chosen to be U£l = 17.2 MeV. Only the states dominated by the QfiVt configuration 
are presented. See the text for the discussion of other states. The experimental values are taken 
from Refs. p5p5|. 



6 He 


k Q = l 


k Q = 0.6 


Exp 


Eb 


33.230 


27.953 


29.269 


Ad 


1.629 


1.707 


1.72 




2.196 


2.317 


2.59 


eM^) 











eM^) 


2.577 


2.449 


1.80 


eM^) 


6.339 


5.638 


N/A 


eM,V 


7.140 


6.324 


N/A 


eMA) 


9.565 


8.355 


N/A 



TABLE V. Experimental and calculated binding energy, point proton and neutron radius in 
fm, and excitation energies E X (J W ,T) for 6 He. All energies are in MeV. The presented results are 
for 6h£l and are performed as described in Table [V. The experimental values are taken from Refs. 
|2|[27|. 
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